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Análise linear elástica pelo método dos elementos 
finitos do provete DCB usado em mecânica da fractura 


SUMÁRIO 


A deformação linear elástica de um provete DCB 
(double cantilever beam) foi analisada usando o mé- 
todo dos elementos finitos. Os resultados foram com- 
parados com medições experimentais levadas a cabo 
em provetes de alumínio. A relação entre a abertura 
da fissura à face do provete, medida com um transdu- 
tor do tipo COD clip gauge, e o deslocamento nos 
pontos de aplicação da carga, relevante para um con- 
junto de situações referidas, é apresentada. 


1. INTRODUÇÃO 


Os provetes mais comuns em estudos experimen- 
tais de Mecânica da Fractura, tais como o de flexão 
em três pontos ou o CT (compact tension) apresentam, 
a carga constante, um aumento do factor de intensi- 
dade de tensão K com o comprimento da fissura a. 
Porém, se um provete for concebido de tal maneira 
que a sua compliance (deslocamento dos pontos de 
aplicação da carga/carga) seja uma função linear do 
comprimento da fissura, então o factor de intensidade 
de tensão para uma dada carga será constante e inde- 
pendente do comprimento da fissura, como resulta da 
seguinte relação da Mecânica da Fractura 


K=P| 28 101=05 da tm) 

na qual C é a compliance, E o módulo de Young, v o 
coeficiente de Poisson e B a largura da fissura. 

Um tal provete foi introduzido por Mostovoy et 

al [11], que salientaram que para uma geometria deste 


P. M. S. TAVARES DE CASTRO (') 
|. C. RADON (2?) 
L. E. CULVER() 


ABSTRACT 


The Finite Element Method was used in the linear 
elastic analysis of a DCB (double cantilever beam) 
specimen. The present results are compared with expe- 
rimental measurements carried out in aluminium speci- 
mens. 

A relationship is established between the crack 
opening displacements, which have been measured 
with a COD clip gauge, and the displacement of the 
points where the load is applied. This relationship is 
established for a relevant range of crack lengths. 


tipo, consistindo basicamente num par de vigas, obtém- 
-se apreciáveis deslocamentos a cargas relativamente 
baixas. Consequentemente são facilitadas medições ri- 
gorosas da compliance e de dC/da. 

De acordo com a discussão apresentada em [1], 
a seguinte relação 


(3a</hº) + (1/h) = const (2) 


na qual h representa a altura do provete na extremi- 
dade da fissura define o contorno que satisfaz a con 
dição dC/da=const. Uma descrição completa deste 
provete é oferecida por Radon et al [21]. 

O provete DCB tem sido usado muito frequente- 
mente para diferentes estudos de Mecânica da Frac- 
tura, tais como determinação da resiliência em estado 
de deformação plano [3], [4], correlações de veloci- 
dade de propagação de fissuras de fadiga com para- 
metros da mecânica da fractura linear elástica [3 a 6], 
estudos de propagação de fissuras de fadiga em con- 


() Assistente, Faculdade de Engenharia do Porto, Depart. de Eng. Mecânica. Trabalho realizado no Imperial College, Londres 


Endereço actual: Cranfield Inst. of Technology, Depart. of Materials, 


(7 Research Fellow, Imperial College, Londres 
(7 Senior Lecturer, Imperial College, Londres. 


Comunicação apresentada ao Congresso da Ordem dos Engenheiros, Lisboa, Novembro 1977. 
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dições elasto-plásticas [7], estudos de corrosão [8], 
bem como estudos de propagação de fissuras em con- 
dições de fluência [9], [10]. 

Em alguns destes estudos foram feitas tentativas 
de obtenção de um modelo analítico da deformação 
ds provete. Uma destas tentativas consiste em tratar 
os braços do provete como vigas estudadas pela teoria 
das vigas da Resistência de Materiais [10], [111]. 
Isto é, no entanto, criticado por alguns autores [12], 
[13] que defendem que tal análise pode conduzir a 
deformações substancialmente diferentes das observa- 
das experimentalmente. Um engenhoso refinamento 
deste tipo de análise é descrito em [11], onde o autor 
acrescenta às flechas calculadas pela teoria das vigas 
uma parcela que toma em consideração o efeito de 
singularidade na extremidade da fissura, calculada su- 
pondo que a viga poderia rodar em torno da extremi- 
dade da fissura, estando encastrada na extremidade do 
provete. Foram obtidos resultados satisfatórios, mas 
esta achega carece de justificação teórica. 


2. PRESENTE ANÁLISE 


Dado o contorno complicado do provete, uma aná- 
lise pelo método dos elementos finitos parece ser par- 
ticularmente aplicável; isto foi sugerido em [13], e 
levado a cabo em [14], onde os autores analisaram a 
influência da geometria da extremidade da fenda na 
compliance do provete, usando os elementos TRIM6 
do sistema ASKA. 


Foi decidido desenvolver um programa de ele- 
mentos finitos para analisar uma geometria DCEB, 
Fig. 1, correntemente utilizada no Imperial College, 
Londres, tratando como primeiro passo o comporta- 
mento linear elástico. 


É frequentemente necessário determinar a com- 
pliance de um provete com base em medições que 
não são realizadas directamente na linha de aplicação 
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da carga, tal como sucede em estudos de corrosão [13]. 
Outra situação do mesmo género pode resultar do 
seguinte: as determinações de compliance são normal- 
mente levadas a cabo em máquinas de ensaio equipa- 
das com registadores de carga e do movimento do 
travessão, como é relatado, por exemplo, em [9], [11]. 
[13] em que foram utilizadas máquinas Instron; é de 
esperar, porém, que medições do deslocamento dos 
pontos de aplicação da carga realizadas com transdu- 
tores aplicados no provete sejam mais rigorosas. Na 
presente geometria, Fig. 1, surge a dificuldade de po- 
sicionar o transdutor (clip gauge) na linha de aplicação 
da carga. Portanto, este tem de ser colocado em apoios 
convenientes (knife edges) afastados daquela linha. Foi 
então decidido analisar pelo método dos elementos 
finitos a relação entre o deslocamento medido pelo 
transdutor, e o deslocamento que realmente interessa 
para a determinação da compliance, ao longo da linha 
de aplicação da carga. 


3. O PROGRAMA UTILIZADO 


O método dos elementos finitos é hoje uma ferra- 
menta comum para a solução de problemas de mecã- 
nica dos meios contínuos, bem documentada em vá- 
rios livros de texto de diversa extensão, como por 
exemplo [15], [16]. 

À presente análise foi levada a cabo com base 
num programa de elementos finitos disponível no Im- 
perial College, o FLAP3 introduzido por Slater [17]. 
Resumidamente, o programa analisa elasticamente os 
pequenos deslocamentos num corpo plano e isotrópico 
da espessura uniforme, sob condições de estado de 
tensão ou deformação plano. O plano do corpo é divi- 
dido em triângulos, que constituem elementos de es- 
pessura unitária, e o programa sintetiza os elementos 
compondo o corpo de tal modo que existe equilíbrio 
e compatibilidade de deslocamentos nos nós dos ele- 
mentos. 


O O e + 
= = = : 
H 


FIG. 1 


Provete DCB utilizado no Imperial College, Não se representam entalhes no plano de propagação da fissura (plano de 
simetria) 
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As características de deformação dos elementos 
são tais que o deslocamento varia linearmente em cada 
elemento, o que assegura que eles se ajustarão perfei- 
tamente quando deformados. As tensões são calculadas 
a partir dos deslocamentos nodais, e se u e v repre- 
sentam deslocamentos que variam linearmente no ele- 
mento nas direcções cartesianas x e y, então as compo- 
nentes da deformação 


dv 
dx 


dv du 


(3) 


dy dy 
são constantes em qualquer elemento. As tensões q,, 9 


T,.. São dependentes de e,, e, e, como segue: 


x 


v* 


Xv 


O, ; 1 y 0) e, 
E À ss] 9 1 0 e, 4 (4) 
: E 0 O (1-v)/2 e 


e serão igualmente constantes em qualquer elemento. 

O método de solução da equação de rigidez, no 
presente programa, é a iteração de Gauss-Seidel, com 
sobre-relaxação. 


4. GERAÇÃO DA REDE 


A rede foi concebida por forma a definir o con- 
torno do provete o mais fielmente possível, o que é 
deveras simples quando são usados elementos trian- 
gulares, visto que o contorno pode ser aproximado por 
uma série de pequenos segmentos de recta, que cons- 
tituem os lados dos triângulos. Outras considerações 
têm porém de ser levadas em conta, como a necessi- 
dade de ter elementos tão equiláteros quanto possível, 
e a vantagem de um método automático de numerar 
elementos. 

A rede gerada é apresentada na Fig. 2. Os seguin- 
tes pontos devem ser sublinhados: 
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Técnica de numeração dos elementos 


A fissura é simulada impondo dois diferentes 
conjuntos de condições fronteira ao longo do 
eixo dos x: os nós dos elementos correspon- 
dentes ao comprimento da fissura em estudo 
têm liberdade de deslocamento em ambas as 
direcções x e y, excepto o nó correspondente 
à extremidade da fissura; este, juntamente 
com os outros nós do plano de simetria per- 
tencentes à parte não fracturada do provete, 
apenas se podem deslocar na direcção x. 

Como estaremos tratando deslocamentos me- 
didos nos apoios do transdutor (knife edges), 
Fig. 4a, o provete é suposto ter o aspecto 
representado na Fig. 4b, o que permite a 


-—- a o ál 
1483 1607 1348 


Provete DCB. Rede de elementos finitos 
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1. Apenas metade do provete foi considerada, 


tendo em vista a sua simetria, 

2. Com excepção de 2 dos 512 elementos que 
constituem a rede (elem. 1 e 12, Fig. 3), esta 
é uma tira constituída por linhas de dez ele- 
mentos, enrolada em torno do orifício do 
pino de aplicação da carga. Naturalmente, 
esta concepção permite um muito expedito 
processo de numeração, Fig. 3. 
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comparação entre a análise pelo método dos 
elementos finitos do deslocamento do ponto 
A, Fig. 4b, com o deslocamento medido ex- 
perimentalmente nos apoios (knife edges) pelo 
transdutor. Esta aproximação justifica-se invo- 
cando que numa viga cantilever carregada 
num ponto situado entre a secção de encas- 
tramento e a extremidade livre, o comprimen- 
to livre da viga não está sujeito a qualquer 


ia a) 


[b] 


FIG. 4 


Provete [a] e idealização para análise pelo método 
dos elementos finitos (b) 


tensão, e deflecte linearmente segundo a tan- 
gente à deformada no ponto de aplicação da 
carga (desprezando o efeito do peso próprio). 


As coordenadas da rede foram obtidas a partir 
d: um desenho à escalada 2/1, usando um CADMAC 
digitizer. Esta máquina lã directamente coordenadas 
x, Y. que imprime em cartão ou fita perfurada. 


5. RESULTADOS E DISCUSSÃO 


O programa de elementos finitos foi corrido num 
computador CDC 6400. 

Antes de tratar o problema em análise, o programa 
foi testado numa situação típica, nomeadamente a de 
uma viga cantilever de secção rectangular, carregada 


a; 


| 
PáirdiVá 
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Primeiro, conforme sugerido por Gross, Roberts 
e Srawley [18] o programa foi corrido na versão es- 
tado de tensão plano, visto que no provete apenas a 
região próxima da extremidade da fissura se aproxima 
do estado de deformação plano, enquanto o restante 
material se encontra sujeito ao estado de tensão plano. 

Em segundo lugar, a carga vertical aplicada foi 
distribuída pelos nós situados no semicirculo superior 
do orifício do pino de aplicação da carga. 

A base de comparação para a análise pelo método 
dos elementos finitos do provete DCB consistiu em 
medições efectuadas com um transdutor (COD clip 
gauge) aplicado em apoios apropriados (knife edges, 
Fig. 4a). Estes apoios foram colados ao provete usando 
ARALDITE RAPID, com uma montagem conveniente 
para os fixar na posição correcta até à secagem da 
cola. O transdutor foi calibrado usando um calibrador 
INSTRON com uma resolução de 0.0025 mm, e O 
respectivo sinal foi registado num registador BRYANS. 

As medições foram obtidas em provetes de alu- 
mínio RR58, durante o programa de ensaios de propa- 
geção de fissuras de fadiga descrito em [41]. 

Quanto ao programa do método dos elementos 
finitos, este foi corrido para diferentes comprimentos 
da fissura, usando as propriedades mecânicas do alu- 
mínio RR5B, e os deslocamentos obtidos foram corri- 
gidos para tomar em consideração a espessura dos 
provetes, 0.95 cm. Na Fig. 6 são comparadas medidas 
do transdutor e resultados do método dos elementos 
finitos. Foi obtido um acordo satisfatório, apresentan- 
do-se os valores experimentais ligeiramente mais ele- 
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FIG. 5 


Viga cantilever. Rede de elementos finitos 


verticalmente na sua extremidade, Fig. 5. Foi encon- 
trada uma diferença de 8 % entre as flechas calculadas 
por aplicação da teoria das vigas da Resistência de 
Materiais e as flechas calculadas usando o programa 
de elementos finitos, a primeira análise apresentando 
os valores mais elevados. Além disto, foi verificado 
que a razão entre as flechas de dois pontos fixos, 
obtida por aplicação da teoria das vigas, ou pelo mé- 
todo dos elementos finitos, diferia menos de 1 %. 

Quanto à análise do provete DCB, dois pontos 
preliminares devem ser referidos: 
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vados que os valores resultantes da análise de compu- 
tador. Atribuindo maior confiança às medições experi- 
mentais, a diferença obtida pode ter a seguinte justifi- 
cação: conforme Zienckiewicz [151], é possível mos- 
trar que na ausência de quaisquer tensões ou deforma- 
çoes iniciais, como no caso presente, o método dos 
elementos finitos subestima o valor da energia de de- 
formação, constituindo uma solução tipo «lower bound». 

A razão do deslocamento na posição do transdutor 
(knife edges) e o deslocamento na linha de aplicação 
da carga, de interesse prático para situações já descri- 
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FIG. 6 


Deslocamento dos apoios do transdutor, ponto A. Medidas expe- 

rimentais e resultados do programa de elementos finitos para a 

mesma carga aplicada (1340N). (Comprimento da fissura, a, medido 
desde a linha de aplicação da carga) 


tas, é apresentada, como função do comprimento da 
fissura, na Fig. 7. Deve ser apreciado que o aspecto 


——— — — + 
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Comprimento da fissura q im | 
FIG. 7 


Análise pelo método dos elementos finitos da razão entre o des- 

locamento medido pelo transdutor (COD clip gauge) e o deslo- 

camento dos pontos de aplicação da carga, em função do com- 

primento da fissura, a. (a. medido desde a linha de aplicação 
da carga) 


qualitativo deste resultado era previsível, dado o com- 
portamento comparável a uma alavanca do provete. 

Os valores numéricos referidos neste relatório es- 
tão coligidos na Tabela |. 


6. CONCLUSÕES 


Uma descrição satisfatória do comportamento de 
um provete DCB foi obtida com uma análise pelo mé- 
todo dos elementos finitos usando elementos CST 
(constant strain triangles). 

Após discussão da importância prática do conhe- 
cimento da razão entre o deslocamento à face do pro- 
vete e o deslocamento nos pontos de aplicação da 
carga, aquela razão é apresentada como função do 
comprimento da fissura. 
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ANÁLISE PELO MÉTODO DOS E. F. DO PROVETE D. C. B. 


N.º de elementos CST: 512 
Carga aplicada: 1340 N 


N.º de nós: 306 
Estado de tensão plano 


E=0.75x10' Nm v=0.33 
a) Output da análise por E. F. 
N.º nós Factor N.º ciclos Va VB Va 
imobil s. relax, Toler iteração Ve 
(em) seg. Y (cm) (em) 
td | 16 1.992 -1000E-4 833 1223E-1 .5320E-1 1.3589 
8.66 14 1.993 2T74E-4 900 9349E-1 TIO3E-1 1.3162 
10.34 12 1.995 .2067E-4 900 1222E 0 9542E-1 1.2807 
12.18 10 1.996 .1866E-4 900 1540E O 1230E O 1.2521 
14.19 8 1.997 .2859E-4 900 A953E O 1591EO0 1.2272 
b) Deslocamento do transdutor em função de a 
MEDIÇÕES EXPERIMENTAIS ANALISE POR E. F. 
a Va a Va (mm) 
(em) (mm) (cm) ' .. 
7.15 1.65 7.14 J23x2x(1/0.95) = 1.520 
9.50 2.49 8.66 935 x 2 x (1/0.95) = 1.969 
10.30 2.60 10.34 1.222 x 2x(1/0.95) = 2.574 
11.35 3.01 12.18 1.540 x 2 x (1/0.95) = 3.230 
12.80 3.50 14.19 1.953 x 2 x (1/0.95) = 4.113 
DO. E E 
E x2, visto que apenas metade do provete foi analisada, 
O X(1/0.95), visto que a espessura do provete original não era a unidade, mas sim 0.95 em, 
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ANO LV 


VOLUME XLI 


Notas sobre o traçado do diagrama psicrométrico do ar 


RESUMO 


Neste trabalho analisam-se os problemas que surgem 
quando se traça um diagrama psicrométrico do ar num 
traçador de gráficos acoplado a um computador. O dia- 
grama que se apresenta é o diagrama entalpia — humi- 
dade especifica com a linha de 0º C na posição vertical. 
Confrontam-se algumas expressões de propriedades a 
representar no diagrama e apresenta-se uma expressão 
para o cálculo do factor térmico ambiente. Apresentam-se 
ainda algumas considerações práticas e úteis para o 
traçado. 


1. INTRODUÇÃO 


O diagrama psicrométrico do ar é um instrumento 
fundamental para o cálculo e dimensionamento das ins- 
talações de ar condicionado, frigorificas, de secagem e 
torres de arrefecimento. Por outro lado é também um 
facto que a maioria dos diagramas disponíveis neste 
momento nem sempre concordam uns com os outros e 
além disso apresentam alguns erros que, embora com 
origem na necessidade de simplificar o traçado, condu- 
zem a não poucas confusões, das quais a mais comum 
é a verificada entre os vários conceitos de temperatura 
húmida. 

O traçado num traçador de gráficos acoplado a um 
computador ou a um calculador permite que, para além 
de uma excelente precisão de marcação dos pontos, seja 
possivel traçar mais correctamente as linhas que aparen- 
tando ser rectas o não são (temperatura húmida e volume 
especifico). Isto deve-se a que estas linhas exigem o 
cálculo dum maior número de pontos e portanto um 
maior volume de cálculo. Este é facilmente suportado 
mesmo por um pequeno calculador. 


2. EQUAÇÕES FUNDAMENTAIS 
2.1. Constantes utilizadas 
Constante universal dos gases perfeitos — R,=8314,3 J 


kmole”! K”! 


Massa molecular do ar seco—m,=28,97 kg kmole”” 


Manuscrito recebido para publicação em 185/77. 
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ABSTRACT 


This paper presents the principal problems occurying 
when drawing a psychrometric chart on a plotter coupled 
to a computing system. In this case the chart to be 
drawn was a distorted enthalpy — moisture content chart 
in order to present dry bulb temperature and moisture 
content as the principal axes and a vertical 0º C dry bulb 
temperature line. A few expressions for the properties 
are compared and a new expression for the room factor 
calculation is presented. A few pratical expressions which 
allow easy plotting of the diagram are discussed. 


Massa molecular do vapor de água—m, = 18,02kg 
kmole”! 


Constante de gás perfeito do ar— R,=287,00) kg! K'! 
Constante de gás perfeito do vapor de água —R, = 


= 461,39] kg! K 
Pressão atmosférica — P,, = 101325 Nm? 


2.2. Volume específico 


Considerando o ar como um gás perfeito 
Pa Va=R, T, 
A pressão parcial do ar seco é a diferença entre 


a pressão total atmosférica e a pressão do vapor de 
água existente na atmosfera e portanto 


BR. (t+273,15) 
ER qi (1) 


2.3. Pressão de vapor de água 


Como é evidente a pressão de saturação do vapor 
de água a uma dada temperatura é função apenas dessa 
temperatura. 
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Para temperaturas Celsius positivas recolheram-se 
duas expressões. Segundo Jones [2] que cita [5] a pres- 
são de saturação do vapor é 


log P,=33,59051 —8,2 log (t+ 273,15) + 0,0024804 x 
(427315) 3142,31 (2) 
o t+ 273,15 | 


expressão válida para o intervalo entre 0ºC e 100º C 
e no qual apresenta o erro relativo máximo de 0,013 % 
a OC em relação aos valores de [3]. 


Por seu lado Bull [1] sugere a expressão 


ia + 2,7857 (3) 
237,3+t.| | 


log Fi = 
que no intervalo de OC a 100ºC apresenta um erro 
relativo máximo de 0,815% a 100"C. No anexo 1 é 
possivel comparar os erros relativos de cada uma das 
expressões propostas e concluir que (2) é bastante pre- 
cisa e (3) é muito razoável especialmente entre 0º C 
e 60 C. Aliás, os desvios padrão calculados em pontos 
equidistantes de 5'C são de respectivamente 2,506 e 
239,488 para (2) e (3) no intervalo entre 0º C e 100º C 
e de 0,399 e 1,969 no intervalo entre 0º C e 60º C. Neste 
último intervalo os erros relativos máximos são para 
(2) e (3) respectivamente 0,013% a 0ºC e 0,048% a 
45 C. No diagrama psicrométrico é muito rara a repre- 
sentação para além de 60" C e portanto não será a pre- 
cisão da expressão que condicionará a escolha mas 
sim a simplicidade e a rapidez de cálculo da função e 
da sua inversa, Perante isto é claro que a expressão 
a aplicar será (3). 
Para temperaturas Celsius negativas a pressão de 
saturação do vapor de água em contacto com gelo é 
segundo Bull [1] 


9,5 t 
— + 27857 (4) 


og Pis 
5 265,5+t 


e segundo Jones [2] citando [4] é 


log P. = 125380977 suo] 5 
o v= tz, e — ee 
q t+ 273,15 63) 


Comparando com os valores tabelados em [3] com 
valores espaçados de 10º C no intervalo de —40ºC a 
0º C obtém-se valores dos erros relativos que se podem 
apreciar sob forma gráfica no anexo 2. Nesse inter- 
valo o erro relativo máximo e o desvio padrão são res- 
pectivamente 0,565% a — 40º C e 0,240 para a expres- 
são (5) e 2,132% a — 40º C e 0,572 para a expressão (4). 
Assim, e por razões idênticas às apontadas acima, ao 
escolher (3) para temperaturas Celsius positivas, esco- 
lheu-se (5) para temperaturas Celsius negativas. 


24. Humidade especifica 


A humidade especifica do ar é a razão entre a massa 
de vapor de água e a massa de ar seco, 
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tw) = MAM, 


Estas massas ocupam o mesmo volume V e sendo 
va & v, OS volumes especificos do ar seco e do vapor 
respectivamente 


V 

| V Va 
V Ve 
W 


e se pudermos considerar o ar seco e o vapor de água 
como gases perfeitos 


RA P, 
= —— —— 
AR; P, 


e como py =P, + Pa, 


R P 
EA E e ie (6) 
à, Par Pv 


A expressão de [6] é inútil pois a correcção que pretende 
introduzir não tem qualquer efeito na posição das linhas 
e complica o cálculo e o traçado. 


2.5. Temperatura de bolbo húmido 


Mais correntemente denominada temperatura húmida 
esta temperatura pode ter normalmente três significados 
distintos: temperatura húmida de saturação adiabática 
(temperatura húmida termodinâmica) e temperatura de 
bolbo húmido protegido e não protegido. 

A temperatura húmida a representar no diagrama é 
a última citada pois é a que se determina normalmente 
com um termômetro cujo bolbo está envolvido numa 
gaze húmida. 

Segundo [1] existe uma relação entre a pressão do 
vapor e a temperatura húmida, que é 


P,=P,—A Pa (t—t) (7) 
O valor da constante A depende da temperatura 


húmida e do termómetro utilizado. Tem os seguintes 
valores 


Tipo a e tt>0 
Não protegido 594x104 KT 6,66xX104K) 
Protegido («screen») 7,20x10* K!  7,99x104 K 


2.6. Entalpia especifica por kg de ar seco 


O estado de referência da escala de entalpia é o 
estado em que o ar se encontra seco (não há qualquer 
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vapor de água misturado) e à temperatura de 0ºC. 
Tomando este estado como estado de referência a ental- 
pia do ar seco a uma temperatura t é 


hac= 1,009 t 


A entalpia do vapor de água saturado pode ser 
aproximada pela expressão 


h,.=2500,8 + 1,84 t 


Se atentarmos num diagrama de Mollier ou nas 
tabelas de sobreaquecimento para o vapor de água, ve- 
mos que, para as temperaturas e pressões de saturação 
que se verificam nas zonas de aplicação do diagrama 
psicrométrico, a entalpia é praticamente só dependente 
da temperatura e tem um valor idêntico ao da entalpia 
de saturação, isto é, a baixas pressões o vapor de água 
comporta-se como um gás perfeito desde a saturação. 

Assim e existindo a humidade específica w, a ental- 
pia específica da massa de ar húmido por kg de ar seco 
será 


h=ha+w h, 
ou seja 


h=t (1,009+1,84 w)+2500,8 w (8) 


2.7. Humidade relativa 


É a razão entre a pressão que o vapor de água 
exerce a uma dada temperatura e a pressão que exer- 
ceria se estivesse saturado a essa mesma temperatura. 
Expressa em percentagem é 


E 
9 =5— X 100% 


vs 


3. O DIAGRAMA NUM PLANO h,tw 


De (8) pode concluir-se que num diagrama h,w, as 
isotérmicas são rectas não paralelas entre si e de coefi- 
ciente angular 


“dh -184t+ 25008 
dty j 


No que diz respeito às isocóricas se em (6) expli- 
citarmos P, e substituirmos em (1) conclui-se que 


Ra 
v , 
t= E Pa = AA 
Edita 
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que introduzido em (8) conduz a 


Va Pu / 1,009+1,84 W 
h = + 1998,204 w — 275,60835 
R, BR, 
+ w 
R, 


que não é a equação de uma recta. A sua derivada é 


— 1,009 


+ 1998,204 


O afastamento da linearidade é causado pelo termo 


Va Pa 1,009+1,840 a 
AR, A 
t) 
qm E 


que se torna de maior peso na expressão quando os 
outros termos se anulam, isto é, quando 


275,60835 Gr Sac 
W=-—————— = 0, 
1998,208 9y a 


nestas circunstâncias a entalpia é igual a 


Quando (9) toma um valor constante por se des- 
prezar tw (w=0), h=h (two, v,) é uma função tal que 


x 1,622 


o que representa um erro de 2,41 % no caso extremo. 


9h 
Se, por exemplo, se calculassem os valores de 
dw [Ya 


para a isocórica de 0,86m' kg obter-se-ia 2063,6 
kilkg, kg, 'J! para w=0 e 2060,6 kJ(kg, kg, —!) ! para 
w=0,015 kg, kg, ! (ar saturado). As isocóricas podem 
então ser traçadas por meio de dois segmentos de recta 
definidos pelos pontos de ar seco, de ar saturado e de 
humidade especifica média destes dois, duma mesma 
isocórica. 


A equação das isotérmicas húmidas obtém-se expli- 
citando P, na expressão (6) e substituindo o seu valor 
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em (7). Explicitando t na equação obtida e substituindo 
em (8) virá 


(1,009 + 


+ 1,84 ) + 2500,8 w 


que pode ainda escrever-se sob a forma 


h=[ 1009] 4 P'y ms: Ra 4 
at A RH, 
q | 
PE. 
+tw| 1,84] E + -— — | + 2500,8 ni 
| AP A 
% f 
| a 
— ——— | 10090 + 184 | Pu 
A, R. : , r R 
E + tas w 
A, ) -— “ A ma 


em que o último termo é o responsável pela não lineari- 
dade da função. No entanto, se o substituirmos pelo valor 
que toma quando w=0, verifica-se que o erro cometido 
cresce quando to aumenta. Esse erro é de cerca de 
1,67% quando w=0,1 kg, kg, ! e o erro cometido ao 
determinar h(to,t) será menor e, só muito excepcional- 
mente, igual a 1,67 %. Se se calculassem os valores de 
| dh 
MET 
kl(kg,kg, ')! para w=0 e 1283 kl(kg,kg, ')! para 
vu = 10,64 10º kg, kg! 

Podemos então proceder da mesma forma que pro- 
cedemos para as isocoras, isto é, aproximar a linha tra- 
cando 2 segmentos de recta, o primeiro entre o ponto 
de saturação e um ponto médio, o segundo entre este 
ponto e o ponto em que a linha de temperatura húmida 
atinge a linha de humidade específica nula. 


: nara t' = 15 €C, teriamos o valor de 139,2 
E 


4. TRAÇADO DE UM DIAGRAMA t,ty 
4.1. Os factores de escala 


O diagrama t tw é, na realidade um diagrama h, tw 
que sofreu uma distorção de modo a que a linha de 
O0ºC de temperatura húmida seja vertical aparentando 
o eixo das abcissas ser o eixo das temperaturas secas. 
Na realidade as linhas de temperatura seca positiva têm 
um declive positivo tanto menor quanto maior a tempe- 
ratura e as de temperatura seca negativa têm declive 
negativo que decresce com a diminuição da tempera- 
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tura. À distorção é obtida pela rotação das linhas de 
entalpia especifica no sentido directo, mantendo-se hori- 
zontais as linhas de humidade especifica. As isentálpicas 
continuam assim a ser paralelas. 

O diagrama tem uma área útil de comprimento x, e 
altura y, tal como a fig. 1 mostra sendo nosso objectivo 
o traçado do diagrama psicrométrico entre as tempera- 
turas t; e t; e as humidades o e twW,.« Como as linhas 
de temperatura seca negativa se inclinam para a esquer 
da à medida que a temperatura baixa, é o ponto de 
saturação a t, que ficará mais à esquerda do diagrama, 
representado pelo ponto 3 da fig. 1 (tsy=t, e w;=to'(t)). 
Hor argumento semelhante é fácil concluir que o ponto 
mais à direita é o de temperatura seca t, e humidade 
especifica tw... OU seja o ponto 5 tal que t;r=t; e 
or = Una 

Considerando o ponto 8 (tg=0C e wWr«=UWs) é 
imediato que a sua entalpia especifica é 


hg = 2500,8 15 
e a entalpia especifica do ponto 5 é 
hs=t, (1,009 + 1,84 t95) + 2500,8 tos 
cuja diferença é 


hs—hg = ts (1 009 + 1 84 ts) 


que na horizontal equivale a x; Ora, quando w=o tem-se 


1,84 
mera 
1.009 


No diagrama os pontos 9 e 10 estarão então dis- 
tanciados de x, que é a distância entre os pontos O e 7. 
Como t,=0"C, a temperatura do ponto 7 será 


hs — hg 
bo— ta =— =4 


“1007 (14 


De modo semelhante, como x, é a distância entre 
3 e 11 que é igual à distância entre 12 e 13. Ter-se-á 


hy — hz 
ta — ty =———— = 
13 12 1,009 
t, (1 009 + 1,84 tds) + 2500,8 Us — 2500,8 ta 
o 1,009 
ou seja 
t t ETR E 
s=— to = — | + — ty 
es | 1,009 ) 
e portanto 


e da De? o 
= A) 
tg 1 1,009 3 


que dista x, do ponto O 
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FIG. 1 


Assim podemos definir um factor de escala f, para 
a linha para a qual w=0 que relacione comprimentos 
com temperaturas e ter-se-á 


Xo 


= — 


ty—ts 


Também é possivel definir um factor de escala para 
os ordenados e que é 


Para que seja possível marcar um ponto do qual se 
conheçam a temperatura e a humidade especifica há 
ainda que contornar o problema derivado do facto das 
linhas de temperatura seca constante não serem verti- 
cais. Os pontos que tenham humidade especifica dife- 
rente de zero terão um deslocamento horizontal em rela- 
ção ao ponto à mesma temperatura e humidade especi- 
fica nula que depende e é proporcional à humidade 
especifica e à temperatura do ponto que se pretende 
representar. Assim 


x= a, tw 


que será positivo ou negativo conforme t for positivo 
ou negativo, o que está de acordo com o que já foi dito. 
As distâncias x; e xs; são 


XxX =—àa;, t, W3 
X4 = à) t> Us 
cuja soma é 
x3 + x, = à (to tos — t; 03) 


e cuja diferença para xy é a distância entre os pontos 
1 e 2. Esta diferença é f, (t;p—t,) e assim 


x + x =xo— f (tp—ty) 
donde 
xo — fi (tp — ty) 
tz Us — ty 3 


d| = 


TÉCNICA 457 


Como é fácil verificar pela figura 2a a diferença 
entre xy e f;(typ—t,), quando t, é positiva, é x, pelo 
que a, é 

xo — fi; (to — ty) 


t> U)s 


d; = 
Observando a figura 2b e por motivos idênticos a 
diferença entre xy e f; (t;— t;) reduz-se a x; e assim a; é 


xo— H (tp — tj) 
Mo DD >> 
t, 693 


Com estes factores de escala o traçado das linhas 
de temperatura seca, humidadde específica e humidade 
relativa constantes é imediato. 


4.2. Traçado das isocóricas 


Para ar seco com um volume específico v, tem-se 
a temperatura 


Va 


— 273,15 


a 


Quando se tratar de um ponto de volume específico 
do ar seco v, situado sobre a linha de saturação do qual 
se pretende determinar a temperatura há que resolver 
em ordem a essa mesma temperatura a equação 


f=(P—P) vi- (t+ 273,15) R, 


cuja solução é rápida por um método iterativo (Newton- 
-Raphson, por exemplo). 

Querendo calcular a humidade especifica de um 
ponto de temperatura t, e volume especifico do ar seco 
Va bastará calcular P, 


BR, (t,+273,15) 
Pe = Pao Cs 


a 


FIG. 2a 


e calcular t) por meio de (6). No caso contrário (ty co- 
nhecido) ter-se-á 


va R, Par 
PP, = ————m—— 
' R+w R, 
e portanto 
Pat Ya 
t= — 273,19 
Rato RB, 


4.3. Traçado das isentálpicas 


As isentálpicas são linhas paralelas que no plano x,y 
têm a equação geral 


y=ax+b 


A linha que na fig. 1 une 8 a 9 é a isentálpica re- 
presentativa de h=2500.8 tw, e assim a; será 


Para calcular o parâmetro b duma dada isentálpica 
bastará então conhecer o ponto em que ela corta a 
linha w=0 kg, kg, !. 


Quando for caso de calcular a humidade especifica 
de um ponto de entalpia h e temperatura t, ter-se-á 


h — 1,009t 
1,84 t+ 2500,8 
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4.4. Traçados das linhas de temperatura húmida cons- 
tante 


O cálculo do ponto sobre a linha de saturação com 
temperatura húmida t' é imediato pelo que não nos 
debruçaremos sobre ele. Para o cálculo de um ponto 
sobre a linha w=0 é evidente que de (7) se tem 


Em casos em que se conheçam t e t, calcular-se-á 
P. através de (7) e introduzindo este valor em (6) obter- 
-Se-g Li, 


4.5. Traçado das linhas de factor térmico ambiente cons- 
tante 


Segundo [1] é possivel calcular o declive que estas 
linhas têm no papel impondo uma diferença de entalpia 
(ver fig. 3) h;—h;=H entre os pontos 1 e 3. Deste 
modo. 

ha — h 


1 
d;=(x,—x) sen £ = TT f; sen B = TE f;, sen B 
e 
d, H sen 8 
“ sen(a+8) 1,009 ' sen (e. -+ [3) 
+ .f 
E t 


a 
', 


FIG. 3 
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(13; — W5,) f 
Por seu lado sen & = deito Ega e A 


dz 
€ 2 1,009 à n sen (a + 3) 
e a sen E=-——. ty 
o que nos leva F H y—4td» sn É 
separando vem 
sen (x +) f, 1 H 


sen « sen B fo 1,009 ww, 


donde tg a = 
RR “a H 
fé 1,009 ww tg B 
e como tg 8B= — a; (ver dedução do declive das isentál- 


picas) e como tw =tw, 


RP EEE E PESO 
9 1 H 


f> 1,009 tw — tw; as 


O factor térmico ambiente é, por definição, 


como ty = t; virá 


ha = tz (1,009 + 1,84 ts) + 2500,8 3 
ha = t; (1,009 + 1,84 0) + 25008 w; 
cuja diferença é 


ha — ha = (133 — (5) (1 84 & + 2500,8) 


que substituída na expressão de r conduz a 


H 1,84 ty + 2500,8 


3 — Us | —r 
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que substituda na expressão de tg « nos permite concluir 
que 
1 
= gg >>> 
dE 1,84 +25008 1º 1 


É 1,009 t—r as 


que é diferente e mais cómoda que a expressão pro- 
posta por [1]. 

Assim, escolhendo uma temperatura ty é possível 
calcular todos os factores térmicos ambientes em relação 
a essa temperatura. 


5. TRAÇADO DO DIAGRAMA 


Com base no anteriormente exposto foi escrito um 
programa para o calculador Helwlett Packard 9830A exis- 
tente no CTAMFUTL. Este programa permite o traçado 
do diagrama psicrométrico do ar entre duas temperaturas 
e entre o valor zero e um valor escolhido da humidade 
especifica. O traçado de cada uma das seis familias de 
linhas do diagrama é executado por seis partes inde- 
pendentes do programa, accionadas separadamente. Além 
disso o traçado do diagrama pode ser feito num inter- 
valo qualquer de temperatura que não inclua o valor 0º C, 
o que possibilita o traçado de diagramas para aplicações 
específicas, como é o caso de diagramas para uso em 
instalações de secagem, armazéns frigoríficos ou túneis 
de congelação. 


A figura 4 apresenta o diagrama mais comum (que 
é o utilizado em climatização) no intervalo de tempera- 
tura—20ºC a + 45 C e até 0,03 kg, kg, !. As figuras 
5 e 6 apresentam diagramas nos quais não é represen- 
tada a linha de temperatura seca de 0º C. Na primeira 
temos um intervalo de temperatura de 15ºC a 60ºC 
com a humidade especifica máxima de 0,035 kg, kg, |. 
Na segunda representa-se o diagrama psicrométrico do 
ar para o intervalo de temperatura de —40"C a —5'C 
e até à humidade especifica de 0,002 kg, kg, !. 


No canto superior esquerdo dos diagramas encon- 
tra-se o ábaco para os factores térmicos ambientes. 
A cruz que se encontra um pouco para a direita do 
ábaco indica o ponto com o qual se deve unir o ponto 
do ábaco com um valor do factor térmico ambiente 
para se obter o declive das linhas do diagrama com 
esse mesmo factor térmico. Nos diagramas que se 
apresentam as temperaturas de referência para a cons- 
trução dos ábacos foram respectivamente 20º C, 50ºC 
e —20'C. 
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Disjuntores Brown Boveri de 
-r comprimido para instalações ao ar livre, 
para as mais elevadas 


capacidades de ruptura 
e tensões de serviço 


E 
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Instalação de disjuntores ao ar livre, de 245 kV, em 
serviço na Electricité de France, com disjuntores de ar 
comprimido, tipo DLF 245, capacidade de ruptura 12 GVA 


Fornecemos: 
Disjuntores de ar comprimido para instalações ao ar 
livre, para tensões de serviço de 72.5 até 765 kV 
Correntes nominais de 2000 até 4000 A 
Correntes de ruptura até 72 kA 

— Elevado grau de segurança 

— (Construção robusta 

— Amplas possibilidades de adaptação ao 
funcionamento com potências mais elevadas de 
ruptura devido ao tipo de construção por blocos 

— Tempos minimos de disparo 

— Aptos a reengate automático rápido 


Aproximadamente 20 000 disjuntores Brown Boveri de 
ar comprimido instalados ao ar livre garantem 
em 45 países o funcionamento de rêdes de alta tensão. 


Representada em Portugal por: 
Sociedade de Electricidade 


Brown Boveri Limitada 


Rua de Sá da Bandeira, 481-2 , Porto 
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SIMBOLOGIA 


A -— constante da equação das isotérmicas húmidas KT 

a, — factor de inclinação das isotérmicas secas mm (kg, kg, !)! “C7! 

a; — declive das isentálpicas 

f, — factor de escala das abcissas mm K 

b — > » . » ordenadas mm (kg, kga ')! 

h — entalpia do ar húmido kJ) kgs! 

has — entalpia do ar seco » 

h, — entalpia do vapor de água kJ kg! 

hos — » » » » » saturado » 

m, — massa molecular média do ar kg kmole”! 

Mm, —  » » do vapor de água » 

P. — pressão parcial do ar seco no ar húmido Nm? 

P. — pressão atmosférica » 

P. — pressão parcial do vapor de água no ar húmido » 

P, — » » » » . » saturado a t » 

R, — constante média de gás perfeito do ar ] kg! KT 

Ro — » de gás perfeito do vapor de água » 

Ro, — constante universal dos gases perfeitos ] kmole”! K”! 

T temperatura absoluta K 

t — temperatura seca do ar Ke 

t' — temperatura húmida do ar ºC 

Va — volume específico do ar seco mº kg 

vv — > » » vapor de água » 

X, — largura útil para o traçado do diagrama mm 

x, — abcissa do ponto i » 

Yo — altura útil para o traçado do diagrama » 

y; — ordenada do ponto i . 

Letras gregas 

74 — ângulo de factor térmico ambiente rad 

[és — humidade relativa % 

to | — humidade absoluta kg, kg! 

to'(t)) — humidade absoluta de saturação a t=t, » 
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and Transport Properties of Fluids — Sl units, 
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October 1964, pags. 268-278. [5] National Engineering Laboratory Steam Tables, 

[2] Jones, W. P. Air Conditloning Engineering, 1964. 
Edward Arnold (Publishers), London, 2nd edi- [6] Porchier, G., Cours de climatisation, Tome |, 
tion, 1973. E. D. 1. P. A. 1973. 
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Erro relativo das expressões da pressão de saturação de vapor de água a temperaturas negativas 
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Estudo retrospectivo de cimento Portland 
de alta resistência inicial com teores de magnésia 
reputados anormais (*) 


SUMÁRIO 


A Memória dá ênfase à inocuidade da presença 
da magnésia no cimento Portland, pelo menos dentro 
dos limites referidos. 

Alude preliminarmente às propriedades essenciais 
do cimento Portland desde o alvor de sua industria- 
lização as quais, vindo-se apromorando, têm originado 
diversificação de características capazes de satisfazer 
exigências sempre crescentes. 

Faz um retrospecto das propriedades de cimento 
Portland comum fabricado no Brasil, no decorrer de 
anos consecutivos, com propriedades perfeitamente nor- 
mais, em que os índices de resistência crescem com 
os teores de MgoO. 

Baseado no éxito industrial é se levado a pro- 
duzir cimento de alta resistência inicial, mesmo em 
presença de altos teores de MgO, com resultados per- 
feitamente satisfatórios. 

Procura dar uma interpretação aos resultados obti- 
dos e cita alguns experimentadores alheios. 

Finalizando faz uma série de recomendações ten- 
dentes à obtenção de cimento Portland de alta resis- 
tência inicial mesmo em presença de elevados teores 
de MgoO. 


GENERALIDADES 


Desde o início de sua manufactura, então, por 
assim dizer, empírica, tem o cimento Portland, neste 
longo período de um século e meio, vindo sendo con- 
tinuamente sujeito a tantos aperfeiçoamentos e pro- 
gressos na técnica de fabricação, não só pelo emprego 
de equipamentos cada vez mais adequados e eficientes 
como também pela utilização de métodos científicos 
de análise e exame, a ponto de o cimento então pro- 
duzido, pela primeira vez por Joseph Aspdin, possuir 
qualidades deveras medíocres, em contraste desenco- 
rajador com as dos cimentos Portland actualmente pro- 
duzidos. Não que deixasse de comportar-se similar- 


M. A. TEIXEIRA DE CASTRO 
Engenheiro consultor 


SYNOPSIS 


The Study emphasizes that magnesia is innoxious 
to Portland cement within the limits referred into. 

The fact, that the main properties of Portland 
cement since its early manufacture days have been in 
a continuous progress, has led to diversify properties 
capable to mest enlarging requirements. 

It is given a retrospective view of ordinary Port- 
land cement produced in Brasil in the course of con- 
secutive years, properties of which are quite normal, 
the strength increasing with the magnes:ia content. 

In view of that success it was conceived to pro- 
duce high early strength Portland cement with unusual 
magnesia content, properties of which are quite satis- 
factory. 

It is looked for an interpretation of the results 
and a few searchers are cited. 

Finally a few factors, that would enable the pro- 
duction of high early strength Portland cement with 
high magnesia content, are pointed out. 


mente aos actuais, isto é, que deixasse de possuir 
fundamentalmente aquelas propriedades que caracteri- 
zam essencialmente o cimento Portland tais como ma- 
nifestadas, na sua aplicação, pela trabalhabilidade, 
plasticidade, indeformabilidade, inalterabilidade, dura- 
bilidade e elasticidade, mas sua hidraulicidade apresen- 
tava Índices muito inferiores aos que subsequentemente 
iam sendo atingidos. 

De início o que se tinha em mira precipuamente 
era o endurecimento em presença de água, isto é, a 
mera obtenção de hidraulicidade da qual derivariam as 
propriedades acima referidas, sem a preocupação de 


(*) Trabalho apresentado ao 2.º Congresso Brasileiro de Engenharia Quimica, São Paulo, 4 a 7 Julho, 1977. 


Manuscrito recebido para publicação em 13/9/77. 
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lhes serem atribuídos valores. Porém, à medida que 
as obras de engenharia iam tomando vulto, se inten- 
sificando e diversificando, aquelas propriedades, em 
face de exigências consequentemente cada vez maiores, 
tiveram de ser forçosamente levadas a um gradativo 
aperfeiçoamento. Mais tarde, com o advento do con- 
creto armado e ulteriormente do concreto protendido, 
foi tomando incremento a necessidade de assegurar 
resistências cada vez mais elevadas. 

Como decorrência de o cimento Portland ter as 
aplicações mais diversificadas e em meios diferentes, 
seu comportamento está sujeito a variar em presença 
das solicitações a que é submetido e do meio que o 
circunda. Assim, embora fundamentalmente o cimento 
Portland se caracterize por certas propriedades gerais 
maxime a hidraulicidade da qual decorre o endureci- 
mento em presença de água, outras surgem e se ma- 
nifestam, cujos graus de intensidade variam não só com 
o valor intrinseco do cimento mas também em função 
do meio em que é aplicado. 

Dest'arte tiveram pois de surgir exigências espe- 
cíficas conforme os casos a enfrentar nas aplicações 
diversificadas do cimento Portland tais como as que 
se relacionam com a agressividade das águas salinas, 
com as massas ciclópicas de concreto, com as cons- 
truções aceleradas, com as estruturas altamente soli- 
citadas, com a imunidade a certos agregados, etc. 
Quer dizer; tiveram que ir-se imprimindo ao cimento 
Portland, avivando-as, propriedades tais que, sem sair 
da linha mestra derivada da hidraulicidade que lhe é 
peculiar, lhe assegurassem certas características suple- 
mentares que, longe de perturbar o mecanismo funda- 
mental das reacções hidráulicas, antes o tornassem 
mais apto à obtenção das qualidades desejadas na sua 
finalidade específica. Surgem assim, pois, como diver- 
sificações do cimento Portland, além do comum, o 


de baixo calor de hidratação, 

de endurecimento rápido, 

de alta resistência inicial, 

imune às águas salinas, 

imune à reactividade de agregados, 


entre os principais. 


Para enfrentar tal diversificação e consequentes 
exigências o produto terá de ser fabricado em condi- 
ções tais que lhe sejam asseguradas as características 
almejadas, peculiares a cada variedade, produzida, bem 
entendido, uma de cada vez, Como o processo manu- 
factureiro é contínuo, a mudança duma variedade para 
outra numa mesma fábrica pode resultar obviamente 
complicada se se pretender executá-la abruptamente. 

Dependendo fundamentalmente a qualidade do 
produto acabado da sua constituição quimica e estru- 
tural, a instalação existente disponível, com seu apa- 
relhamento mecânico, pneumático e térmico, desem- 
penha todavia papel também preponderante na sua 
consecução. De modo que, tendo de jogar com a ma- 
téria-prima antes do que com qualquer outro factor 


(*) Mátodo de Ensalo Padrão Brasileiro. 
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relativo aos processos inerentes ao seguimento da 
fabricação, não podendo mudar de instalação dum 
momento para outro e bem assim tendo de aceitar as 
matérias-primas tais como ocorrem nas jazidas abas- 
tacedoras, adviria ao fabricante um problema árduo ao 
ter de enfrentar uma mudança brusca e radical na 
qualidade do produto. Assim, uma vez implantada a 
instalação, procura-se no processo manufactureiro che- 
gar, em presença das condições existentes, ao melhor 
desiderato possivel não podendo, sem sacrifício de 
maior ou menor mpnta, ir variando substancialmente 
as caracteristicas do produto a cuja uniformidade busca 
tanto quanto possível cingir-se. Por esse motivo, pelo 
menos em parte, e por não ter havido procura sufi- 
ciegntemente Justificativa, não têm surgido senão em 
escala reduzida no Brasil a fabricação dessas qualida- 
des diversificadas embora a indústria de cimento Port- 
land já aqui perdure há meio século. 

Entre as variedades que têm merecido  parti- 
cular atenção destaca-se o cimento Portland de alta 
resistência inicial, cuja fabricação teve início há umas 
quatro décadas e em torno do qual têm sido efectua- 
dos estudos [71] e surgido discussões. 

O presente trabalho não é senão uma continuição 
cu, melhor, condensação de algo que o autor tem feito 
a esse respeito, baseado em cimentos industrialmente 
produzidos. 


RETROSPECTO DE CIMENTOS PORTLAND COMUNS 
COM TEORES CRESCENTES DE MgO 


Como atrás dito, procura-se em geral, com todas 
as variáveis que entram em jogo no decorrer da fabri- 
cação, chegar a um desiderato considerado óptimo. 
Consegque-se assim, nessa ordem de ideias, melhorar 
gradativa e consideravelmente as propriedades mecã- 
nicas dum cimento Portland comum tais como assina- 
ladas na Tabela | e representadas nos Gráficos 1 e 2. 
Os cimentos designados por A, B, C, D e E corres- 
pondem aos produzidos em anos consecutivos a partir 
do início da fabricação. A ordem alfabética equivale 
à cronológica e os dados abrangem de per si todos os 
obtidos, nos respectivos anos, sobre cimento periodi- 
camente produzido e ensaiado. 

Os métodos de ensaio, sendo anteriores à elabo- 
ração do BM-1 (*), divergem deste mas, para efeito 
de comparação, não afectam em nada as ilações loca- 
lizadas no presente estudo porquanto os ensaios foram 
todos executados em idênticas condições e pelo mesmo 
operador. 

Os ensaios à tração, ao tempo quase universal- 
mente praticados, recairam sobre argamassa destituída 
da plasticidade, de traço 1:3 em peso, sendo a areia 
de rio graduada entre 0,6 e 0,85 mm e a água empre- 
guda ao redor de 10% sobre o peso do material só- 
lido total. Moldagem feita manualmente mediante es- 
pádula em formas metálicas com o feitio de oito, cuja 
secção intermédia, onde se dava a ruptura, media exac- 
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tamente uma polegada quadrada. Esse método, embora 
tanto divergente, se assemelhava ao método adoptado 
pela ASTM nos E. U. A. sob a designação de C 190 [2]. 


Os ensaios à compressão recaíam sobre argamassa 
plástica, de traço 1:5 em peso, cuja consistência, obtida 
numa mesa semelhante à actualmente padronizada no 
MB-1, correspondia a um aumento duas vezes o diã- 
metro inicial da argamassa após 30 pancadas suces- 


sivas, enquanto no MB-1 o diâmetro da argamassa, 
após as pancadas, deve sofrer um aumento de ta 
vezes, quer dizer, a argamassa então ensaiada era mais 
húmida do que a actual argamassa normal. Evidente- 
mente, tratando-se dum traço mais fraco e de consistên- 
cia mais aquosa, o factor água/cimento teria de ser con- 


sideravelmente maior do que o geralmente obtido quan- 
do da adopção do MB-1. A areia empregada era com- 
posta, graduada, com módulo de finura 3,75 e tamanho 
máximo de 9,5 mm. Os corpos de prova eram cubos de 
7,5 cm de aresta. 

Passando-se ao exame dos Gráficos 3 e 4 depa- 
ra-se flagrante ascensão nas resistências mecânicas em 
presença de teores paralelamente crescentes de MgO, 
que variam de 4,3 a 7%, conforme elucidado pela 
Tabela Il, que resume as características químicas dis- 
poníveis relativas aos cimentos a que tais Gráficos se 
referem. Do exame das Tabelas | e Il se infere que 
as resistências aumentam invariavelmente em presença 
de teores crescentes de MgO, sendo essa mesma ten- 
dência acompanhada pelo silicato tricálcico ou pela 


TABELA | 


Características físicas e mecânicas de cimentos Portland comuns constantes da Tabela Il 


DESIGNAÇÃO A B Cc D t 
FIN SI DD rms es UA A 85,3 85,2 86,0 85,5 87,1 
Po Se O a oi a SS Si 1.30 2.10 2.45 2.30 2.20 | 
CONSISTÊNCIA NORMAL, água em peso, % 25,1 24,9 25,3 25,3 25,3 
EXPANSÃO, Le Chatelier, mm ... co ca 3 7 5 4 6 
RESISTÊNCIA, k/cm? | 
TRAÇA é ar som xavier coxa ové 1d — — senso — 15,6 
ES MA DEDO «cas mas ami cu BA — — toi — 19,6 
Areia entre 0,6 e 0,85 mm ... ... 3d — 13,5 16,0 17,5 24,2 
7d 13,4 16,0 18,3 20,7 26,0 
Água, 10 % sobre material seco . 28d 17,4 18,6 20,2 23,1 27,5 
3m 19,1 19,6 20,9 23,6 26,4 
COMPRESSÃO: . sus mas séricos aii cs ZA -— — 75 94 101 
TD OM DURO. a umas cuia a ne 3d -— — 104 122 131 
Areia graduada, mod. 3,75, tama- 
nho máx. 9,5mm .... cc. 7d — — 149 167 181 
28 d — — 223 240 249 
am oe a 282 292 292 
Consistência, alargamento X2 ... 6m — — 295 308 —— 
AQUECE. ass mas ctids Caéi o GRE CaÁ — —. 0,799 0,783 0,786 
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TABELA II 


Características químicas de cimentos Portland comuns constantes da Tabela | 


DESIGNAÇÃO 


ANÁLISE QUÍMICA 


Res. Insolúvel . ... nd. 
Sioz 22,55 
AIZo3 6,77 
FOMOS am uva 2,63 
CARO: sum iqaeço pues 60,92 
MgO 4,31 
SOS: a amerane us pues CRM à 1,51 
Parda ao fogo . ... ... .. 1,38 
RD ços caos ao 100,07 
CONSTITUIÇÃO POTENCIAL (*) 
calculada segundo Bogue 
CAS uam tea 23,1 
ri 47,2 
C3A 13,5 
C4AF 8,0 
CaO.503 2,6 
c3s/C2S “049 


(*) No presente Estudo adoptaram-se fórmulas quimicas simplificadas dos constituintes do cimento Portland, tais como: 


Ca3s — 3Ca0.S102 

C25 — 2C80.5102 

CIA — 3CaO .A!I2O3 

CSAF — 4Ca0 AIZOI.Fe203 


relação C3S/C2S. A única discrepância é uma rever- 
sibilidade que ocorre aos 3 meses com o cimento E, 
facto esse aliás comum aos ensaios à tração em idades 
alongadas. 

Ressalta imediatamente que, pelo menos, a magné- 
sia não constitui um óbice à elevação da resistência 
nem, tão pouco, se afigura um factor adverso que afecte 
a hidraulicidade e integridade do cimento. Com a ele- 
vação de seu conteúdo conseguiu-se melhorar a qua- 
lidade sem que isso Implicasse em modificar subs- 
tancialmente qualquer detalhe no curso da fabricação. 
Apenas se recorreu a um proporcionamento cada vez 
mais adequado das matérias-primas de modo a permi- 
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tir a formação de silicato tricálcico no mais alto grau 
compatível com as condições existentes, tornando-se 
sua elevação tanto mais necessária quanto mais alto o 
teor de magnésia presente, dentro naturalmente de cer- 
tos limites, mas os resultados bem indicam quão longe 
ainda se pode ir no que diz respeito a substâncias que 
foram consideradas outrora não só inertes mas também 
activamente adulterantes e perigosas. E de tal modo 
que, do exame mais atento dos resultados, bem se 
evidencia que a causa determinante dessa melhoria na 
qualidade do cimento, assinalada pela progressiva as- 
censão na resistência, reside não propriamente na 
magnésia em si mas no constituinte de maior poten- 
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cialidade hidráulica do clínquer, em cuja formação 
ígnea a magnésia, dela participando, deve ter benefi- 
camente intervindo. 

Quanto à constância de volume foram os cimen- 
tos naturalmente submetidos ao ensaio então seguido e 
ultasriormente adoptado pela ABNT (*), em 1937, tal 
como consta do MB-1 [3]. Os resultados, embora den- 
tro do limite imposto pela EB-1 [4], afiguraram-se um 
tanto elevados mas releva notar que esse ensaio era 
executado sobre cimentos acabados de produzir para 
os quais, após sazonados ou arejados num curto pe- 
ríodo de tempo, esses resultados praticamente se nuli- 
ficavam. É bem de ver que a causa determinante desse 
fenómeno residia na cal que restara livre em quantidade 
maior ou menor, já que as reacções e combinações no 
decorrer da clinquerização não eram levadas a um 
grau de complementação que mais tarde se tornou 
usualmente praticável pela obtenção de misturas mais 
finamente subdivididas, mais intimamente homogeneiza- 
das e mais intensamente calcinadas. 


CIMENTOS PORTLAND DE ALTA RESISTÊNCIA INI- 
CIAL, COM TEORES ELEVADOS DE MgO 


Em face dos resultados atrás apresentados evi- 
denciava-se a possibilidade de prosseguir na obtenção 
de resistências ainda mais elevadas se bem que não se 
procurasse atingir esse fim simplistamente a partir da 
elevação do óxido de magnésio. Este elemento associa- 
do, em maior ou menor grau, aos calcários, quer no es- 
tado de carbonato sob a forma de dolomita e mais 
raramente de magnésita, quer no de silicatos mais ou 
menos complexos, tem constituído por vezes um en- 
trave na selecção das matérias-primas. Em virtude da 
intensificação e desenvolvimento da indústria de ci- 
mento Portland, vão-se deparando calcários com teores 
cada vez mais elevados daquele elemento de modo que 
a utilização de calcários sucessivamente mais magne- 
sianos vai-se tornando cada vez mais premente. 

Tem-se procurado, como medida de auto-detesa, 
beneficiar certos calcários, principalmente aqueles onde 
a magnésia ocorre sob a forma de silicatos, facilmente 
pulverizados e elimináveis, tendo-se conseguido, me- 
diante britagem, peneiração e lavagem, reduções apre- 
ciáveis que no produto acabado se reflectem em escala 


aumentada da ordem de 1 vezes. Já sob a forma 


carbonatada não é possivel erradicar a magnésia senão 
recorrendo a tratamentos químicos e físico-quimicos que 
redundariam economicamente proibitivos para o supri- 
mento da fabricação de cimento Portland. 

Assim, partindo de matérias-primas usualmente 
empregadas no fabrico de cimentc Portland comum, 
dirige o autor a atenção para a obtenção de cimento 
Portland de alta resistência inicial, com teores de MgO 
invulgares para esta variedade de cimento sem serem 
contudo dispares em relação aos que usualmente se 
deparam no cimento Portland comum. 

O cimento Portland de alta resistência inicial, 
como sua designação mesmo indica, é um cimento Port- 


(*) Associação Brasileira de Normas Técnicas. 


268 


land que deve atingir, nas primeiras idades, resistên- 
cias comparáveis às que se obtêm em idades mais 
avançadas nos cimentos ditos comuns, de tal forma 
que, em geral, o primeiro assegura resistências finais 
num nível ainda mais elevado que o obtido com o 
cimento Portland comum. É, por assim dizer, um ci- 
mento dotado de potencialidade revigorada, isto é, 
levada ao máximo. Para atingir esse desiderato tor- 
na-se pois necessário estabelecer condições de fabrico 
tois que permitam propiciar a obtenção do constituinte 
dotado de mais forte hidraulicidade que se reduz essen- 
cial e esquematicamente, como atrás dito, sob o ponto 
de vista termo-químico, ao silicato tricálcico, cuja pre- 
dominância deve ser inexoravelmente assegurada. É este 
porém um dos constituintes do cimento Portland cuja 
formação se torna tanto mais difícil quanto maior for 
à basicidade da mistura calcário-argilosa que lhe vai 
dar lugar, por realizar-se igneamente sob a forma só- 
lida, sendo não obstante facilitada em concomitância 
re fase líquida que é determinada principal e propor- 
cionalmente pelo óxido de ferro, dentro de certos limi- 
tes e, quiçá, até um certo ponto pela magnésia se bem 
que em menor grau. Conseguiu-se assim com a eleva- 
ção do ferro, a uma temperatura menos alta, redu- 
zindo a quantidade de CaQO livre ao mínimo, obter o 
máximo de silicato tricálcico, tal como apresentado na 
Tabela III. 

Embora residindo essencialmente na estruturação 
mineralógica do clínquer a potencialidade hidráulica 
que ulteriormente se evidencia em contacto com água, 
torna-se necessário demais, para que essa hidraulici- 
dade seja levada ao máximo de eficiência, que às par- 
ticulas em que o cimento se subdivide se assegure um 
máximo de dispersão e reactividade, o que se torna 
praticável pelo maior estado de subdivisão a que o ci- 
mento é levado. 

Assim o cimento Portland, após produzido tal 
como ordinariamente, foi submetido a nova mpagem 
tendo sido necessário, para normalizar a pega, incor- 
porar-lhe gesso adicional. O grau de finura foi levado 
tão longe que os cimentos, tanto os da Tabela Ill como 
os da Tabela IV, atingiram uma superfície específica, 
tal como determinada pelo turbidímetro de Wagner, da 
ordem de 3000 cm*/g enquanto para os cimentos co- 
mumente fabricados a superfície específica era da 
ordem de 1800 a 2000 em?/g. 


Os cimentos constantes da Tabela Ill, cujas resis- 
tências à compressão estão separadamente representa- 
das no Gráfico 5, são dotados de características físi- 
cas e mecânicas que os colocam inquestionavelmente 
no campo dos cimentos de alta resistência inicial. 
Já os constantes da Tabela |V, representados no Grá- 
fico 6, que também foram produzidos para idêntica 
finalidade, ficaram abaixo dos dois primeiros referidos, 
constantes da Tabela Ill. Para estes, porém, que ultra- 
passaram plenamente os valores de resistência exigi- 
dos para o cimento Portland de alta resistência inicial, 
O constituinte férrico é consideravelmente mais elevado 
e a relação C3S/C2S é também em geral mais alta do 


TÉCNICA 457 


